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Abstract 

This paper investigates the question whether a realistic black hole can be in principal similar to a star, having a large but 
finite redshift at its horizon. If matter spreads throughout the interior of a supermassive black hole with mass M ~ 10 9 Mq , 
it has an average density comparable to air and it may arise from a Bose-Einstein condensate of densely packed H-atoms. 
Within the Relativistic Theory of Gravitation with a positive cosmological constant, a bosonic quantum field describing H 
atoms is coupled to the curvature scalar with dimensionless coupling £. In the Bose-Einstein condensed groundstate an exact, 
self-consistent solution for the metric occurs for a certain large value of f, quadratic in the black hole mass. It is put forward 
that £ is set by proper choice of the background metric as a first step of a renormalization approach, while otherwise the 
non-linearities are small. The black hole has a hair, the binding energy. Fluctuations about the ground state are considered. 
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1. Introduction 

The standard knowledge that black holes (BHs) have 
no hair (i. e., they are determined by mass, charge and 
spin only) has recently been tested versus observational 
data from quasars. Schild et al. carefully examine why 
standard models of quasars with the various structural 
elements fail to be compatible with the observations of 
a luminous ring at the innermost edge of the accretion 
disc of the quasar Q0957+561 A,B. The observations 
imply that this quasar contains an observable mag- 
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netic moment, a "hair", which represents strong evi- 
dence that it does not have an event horizon [1] . On 
a different track, one of us presented a class of exactly 
solvable BHs of with the Schwarzschild on the outside 
having with one hair, namely the binding energy [2]. 

A great deal of observational activity is presently un- 
derway to measure with 10 micro-arcsec resolution the 
Sgr A black hole object at the center of our Galaxy, 
because it is widely claimed that detecting the shadow 
of the black hole would prove the existence of an event 
horizon. However the black hole with hair objects ei- 
ther indicated by the Schild et al (2006) observations 
[1] or by the Nieuwenhuizen (2008) theory [2], would 
also be a compact object with a strong gravity shadow 
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that is virtually indistinguishable from a black hole 
event horizon. Indeed, the outer metrics are close to 
the Schwarzschild metric. The two classes of objects 
are in fact best distinguished from one another by the 
existence of a light cylinder effect originating the dusty 
torus, and the clear demonstration of such light cylin- 
der effects in ordinary quasars by Schild et al (2009) 
favors the magnetic over the standard black hole [3]. 

Quasars are probably supermassive black holes 
which occur at the center of most galaxies. Even 
though the one of our Galaxy weighs only about 4 
million solar masses, the typical weight is a billion 
solars, with the present champion at 17 billion. Since 
the Schwarzschild radius of a black hole scales with 
its mass, rs = 2GM/c 2 , its density decays as 1/M 2 , 
implying that bigger black holes have a lower average 
density. This motivates to consider a supermassive 
black hole, say with M — 10 Mq. Then the average 
mass density is on the order of grams per liter, compa- 
rable to air. The fact that they are enormously heavy 
only because of their size, suggests that for these ob- 
jects the physics may be not too difficult or unfamiliar, 
and it invites to study their internal structure as a 
standard problem in a standard theory. 

These estimates of the average density confront the- 
oretical descriptions. On the basis of the Schwarzschild, 
Kerr and Kerr-Newman metrics, black holes (BHs) are 
described as singular objects with all matter localized 
in the center or, if rotating, on an infinitely thin ring. 
Recent approaches challenge this assumption and con- 
sider matter just spread throughout the interior [4-6] . 

Let us point at the following simple connection. Su- 
permassive BHs occur in the center of many galaxies 
and weigh about Mbh = 0.0012 M b ui g e [7]. Let us as- 
sume that they consist of hydrogen atoms and that 
mass and particle number are related as M = vNuin 
with mjv the nucleon mass and some v < 1. Neglect- 
ing rotation, we may compare the H number density 
within the Schwarzschild radius Rs = 2GM/c 2 , i.e. 
riH = 3N/4nr s with the one of densely packed, non- 
overlapping H-atoms, that is, with the Bohr density 
riB = 3/47r<Zo, with ao = 0.529 Athe Bohr radius. This 
yields a mass c 3 (ag/8G 3 mjv) 1/2 = 8.26 • 1O 7 M (we 
take v = 1 here), which indeed lies in the range of 
observed supermassive black holes. The correspond- 
ing mass density is 2.7 kg/liter, but scaling as 1/M 2 
for larger black holes, it is as low as 0.6 g/liter for 
the 17 • 1O 9 M black hole, already below the one of 



air. We shall therefore take as characteristic mass scale 
10 9 M Q = 1.98 • 10 39 kg and write 

M = M 9 -10 9 M Q . (1) 

The assumption of spread-out matter poses the ques- 
tion how matter can have a pressure that allows such a 
state. It was proposed originally by Sacharov that the 
vacuum equation of state p — —p could describe mat- 
ter at superhigh densities [8]. Laughlin and cowork- 
ers assume that matter near the horizon could be in 
its Bose-Einstein condensed (BEC) phase, modeled by 
the vacuum equation of state [4] . Dymnikova consid- 
ers BHs obeying it in the interior, which, however, have 
one or two horizons [5]. Mazur and Mottola take the 
BEC idea over to the interior, and investigate a "gravo- 
star" , of which the interior obeys the vacuum equation 
of state, and which is surrounded by a thin shell of nor- 
mal matter having the stiff equation of state p — +p. 
This solution is regular everywhere [6]. A related sub- 
ject is the description of boson stars [9,10]. 

One of us considered a modification of general rel- 
ativity (GR), the Relativistic Theory of Gravitation 
(RTG) that differs from GR only at very high redshift, 
in particular near the horizon of black holes [11]. This 
modification appears to describe the BH interior like a 
star, be it with strong but finite redshifts. 

We shall study a supermassive BH that exist as a self- 
gravitating hydrogen cloud, in a Bose-Einstein con- 
densed phase. Hereto we employ the Relativistic The- 
ory of Gravitation (RTG), which reproduces all weak 
gravitational effects in the solar system [12,13] as well 
as the ACDM cosmology [11]. We shall extend the re- 
cent approach to this line of research started by one 
of us [2,14,15], where an exact interior solution to the 
supermassive black hole problem was found, in a the- 
ory which is a modification of General Relativity. This 
Schwarzschild-type black hole has a hair, namely the 
binding energy of the matter out of which it is com- 
posed. It has a very large but finite redshift at the hori- 
zon, so there is no sharp horizon, no Hawking radiation 
and no role for Bekenstein-Hawking entropy. From a 
principle point of view, this approach, when extended 
to the rotating situation with an accretion disk, may 
describe the above discussed quasar observations. 

Section 2 discusses aspects of RTG and its general- 
ization of the Schwarzschild metric, together with its 
deformation near the horizon. Section 3 introduces the 
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quantum field theory for the H-atoms and their Bose- 
Einstein condensation. Section 4 presents an exact so- 
lution of the interior metric. Section 5 deals with gen- 
eral solution near the horizon. Section 6 discusses exci- 
tation about the groundstate and the paper closes with 
a discussion in section 7. 



2. Relativistic Theory of Gravitation 

We consider a static metric ds 2 = g^dx^dx" , where 
x^ = (ct,r,8,(f>), with p = 0,1,2,3. Here and in the 
sequel, summation over repeated indices is implied. In 
case of spherical symmetry it has the form 

ds 2 = U(r)c 2 dt 2 - V(r)dr 2 - W 2 (r)dn 2 (2) 

with dO 2 = dd 2 + s'm 2 8d(f> 2 . The metric tensor can be 
read off, = diag((7, -V, -W 2 , - W 2 sin 2 0), while 
g 1 *" — {g^ 1 )^- The gravitational energy density arises 
from the Landau- Lifshitz pseudo-tensor [16], gener- 
alized to become a tensor in Minkowski space-time 
[17,11]. For the metric (2) it takes the form 

, MT 2 / r 2 V'WW 3 , 2 , a 
t = - ^ ; I — \-r V — 5r W 



+ 



87rGr 6 \ 
2r :i VW' 



W 



+ 8rWW' - 2r 2 V - 3W 



(3) 



Let us start with the General Theory of Relativity 
(GTR, GR). The Schwarzschild metric reads in the 
harmonic gauge 

U ° = T = 1 -W = r -lM' W ^r + M, (4) 
Vs Ws r + M 

where M is the gravitational length 



T7 GM 

c 2 



(5) 



The metric is singular at the horizon Ws = 2M, r 
M and it involves the gravitational energy density 



c 4 M d (r + M) 3 (2r + M) 
AnGr 2 dr 2r 3 (r-M) 



(6) 



Its quadratic divergence at M presents an often over- 
looked, physical peculiarity, that induces a negative in- 
finite contribution to the total energy. In other gauges 
the singularity can be moved completely to the origin. 



All by all, the situation is puzzling. For this reason, we 
shall switch to RTG with matter not located at the sin- 
gularity r = 0, but just spread out within the horizon. 

We focus on RTG, which describes gravitation as a 
field in Minkowski space [12,13] and possesses the same 
gravitational energy momentum tensor and thus also 
the gravitational energy density Eq. (3) [2]. It extends 
the Hilbert-Einstein action with the cosmological term 
and a bimetric coupling between the Minkowski (7) 
and Riemann (g) metrics, 



c 4 R 



L = - 167rg ~ PA + |Pbi7Mi'5 M ' y + £mat, 



(7) 



(For phi — it is just a field theoretic descrip- 
tion of GTR.) The dimensionless action is S/h = 
(1/fic) J d 4 x^/—gL, where xo = ct. One has the Ein- 
stein equations 



R" 



1 r> . 
29 K 



8nG 



J tot 1 



1 tot — 1 ' A "•" J bi ' 



(8) 
(9) 



where comes from the matter, from the cos- 
mological term and T£" from the bimetric term. They 
have the dimension of energy /volume and elements 
(p, —pi) of the form 



ptot = p + pa + 



tot 

Pi = Pi - PA + 



Pbi _ Pbi_ 

2(7 2V 
pbi pbi 



pbiT 

W 2 ' 

PbiT 2 

2U 2V ^ W 2 



(10) 



with i = r,0, 4>. The mass density is p/c . It is custu- 
mary to choose the value phi = pa in order to allow a 
Minkowski metric U = V = 1, W — rin the absence of 
matter, because then ptot = p' ot = 0. We shall fix p\ = 
pbi to the observed positive cosmological constant [2]. 
However, historically the opposite choice pA < was 
considered and the cosmological data were described 
by an additional inflaton field [13], so the sign of phi 
is still disputed. We show that the possibility to solve 
a realistic black hole settles that indeed A b i > is the 
physically interesting case. The new point of RTG is 
that goo = U can be very small but still positive. De- 
spite the smallness of pbi, the pu/U term becomes rel- 
evant near the horizon [12,2] and regularizes the sin- 
gularities of the Schwarzschild metric. So we continue 
with p tot w p + Pbi/2U, p- ot w pi + pbi/2U. 
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The pbi term in (7) violates general coordinate in- 
variance and the consistency requirement T£? v = 
imposes the harmonic gauge condition 



u 



V W' 



4rV 

w 2 ■ 



(11) 



The residual gauge group actually would lead to an 
infinite gravitational energy, so its physical subgroup 
is empty, making the solution unique [18] . 

In post-Newtonian approximations and in applica- 
tions to, e. g., stars, the pA and pbi terms are negligible, 
which will bring back the general coordinate invariance 
of GTR. Gravitational radiation in e. g. X-ray bina- 
ries is the same as in GTR. But near the horizon of a 
black hole the large redshift will make the pu/U term 
sizeable, and deeply change the theory. 

The o, i and \ components of the Einstein equations 
read respectively 



1 _ W' 2 

w 2 ~ vw 2 

_L _ w ' 2 

w 2 

u 



2W" V'W' 



vw 2 
w" 

2UV VW 
U'W' 



vw 

U'W' 
UVW 

u' 2 



8tvG 

— Ptot, 



V 2 W c 
8nG tot 

U'V V'W' 



(12) 



+ 



+ 



4U 2 V 4UV 2 2V 2 W 



2UVW 



As usual, the last equation is automatically satisfied by 
energy conservation T!? = 0. The Ricci scalar becomes 

R = ^(-Ptot +Pr 0t + 2 P T) = i?bi + Rm (13) 
R,^^ P f, R m = ^f(-p+Pr+2p^. (14) 

We define the inverse length p\,i by 
p bi = ^|, ,. = ^ = 2.38.10-^,(15) 

with A the cosmological constant. We shall often en- 
counter the combination 

= p ±1 GM = 2 3g 1Q -23 M = 2 3g 10 -i4 M ^ 16) 

C z M0 



2.1. Generalization of the Schwarzschild solution 

The Schwarzschild metric solves the above equations 
in the limit p — > 0, pA = phi — > 0. A more general solu- 
tion within RTG is found as follows. Assume that U, 



V and W are still related in the Schwarzschild manner 
U = 1 - 2M/W, V = W' 2 /U. This solves the two Ein- 
stein equations (12) without matter, p to t = Ptot = 0. 
Inserting this in the harmonic constraint we have 

-f>W- 2 «> + M^r (17) 
Going to the inverse function r(W) this becomes linear, 



W(W - 2M)r" + 2(W - M)r' = 2r 



(18) 



One solution is r = A(W — M). The method of varia- 
tion of constants brings 



A' 



1 

W 



1 



W — M W - 2M 
So with integration constant C one has 



A' = 



-CM 



W(W - M) 2 (W -2M) 
It brings A and from this, choosing A(oo) = 1, 



r = W - M + C 



W 



— M W 
log 

> W -2M 



M 



(19) 



(20) 



(21) 



The horizon U = 0, W = 2M is located at r = oo for 
C > and at — oo for C < 0. Clearly, a small C is 
expected to avoid such peculiarities. 

Near the horizon the bimetric coupling induces a 
deformation of the Schwarzschild metric, which regu- 
larizes its singularity [12,13]. For small values of the 
dimensionless product pbiM, a scaling form was pre- 
sented by one of us [2] , 



■ M 



V- 



1 + ??(e c + C + log y + 2) 
l-r/(eC + C + log77 + 2)' 

w = 



u ■ 

2M 



r/e 



(22) 



r;(l + e<) 2 ' 



l-r/ei - p^iC + wo)' 



Here £ is the running variable and rj ~ PbiM a small 
scale. For pbiM <JC r)e" <C 1 it coincides with the gen- 
eralized Schwarzschild solution (21), with C = —4rj in- 
deed being small. Thus the singularities of (21), U and 
V are smoothly deformed in (22) , due to the bimetric 
term pu/U. 
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3. Quantum field theory in curved space 

Let our H-atoms be described by a scalar, bosonic 
creation field operator 4 



4>(v,t) = V^a;V>i(r)e~ 



iEit/h 



(23) 



where i = {n,£,m}, [a,i,a^] — 8ij and eigenfunctions 
factor as tpi(r) — 4> n e(r)Ye m (8, 4>). The total field [19] 



ipt = i> + i> 



(24) 



is Hermitean. We assume that the two-particle inter- 
action can be replaced by a 5-potential. This leads to 
a quartic Lagrangian density 



L 



full _ 1 f^fl 
mat 2 



A 



■^■(25) 



2fi 2 ' 2"' rl 24nc 1 
where 9 M = d/dx* 1 , (9 M = g^ v d v and A is dimensionless. 
With field dimension [ip] = [v^c/^p], where 




mp 



He 
G' 



(26) 



are the Planck lenght and the Planck mass, respec- 
tively, the Langrangian density has dimension [Lmlt] = 
[hc/£%], i. e., energy density, as it should. 

For a field in curved space the renormalization group 
generates the coupling to the Ricci curvature scalar R 
entering eq. (25) [19]. Its strength £ has to be obtained 
from renormalization arguments. Popular values are 
£ = and £ = ±. In our earlier work [2] we have 
treated it as a phenomenological parameter, here we 
shall argue that its large, mass dependent value is self- 
generated. 

The conjugate momentum field is 



Q r full 



(27) 



The equal-time commutation relation is for the field is 



-93 



(28) 



where n M = S^VU is a time- like unit vector and gz — 
— VW 4 sin 2 8 the determinant of the spatial part of the 
metric. The Hamiltonian density is 



4 In field theory in flat space one usually adds a factor 
1 / yJlEi in the terms of (23); in curved theory space these 
factors are moved to the inner product, see Eq. (78). [19] 



H t = n° t doi>t - lL°I' = (29) 

2 + 2 R)A + 24te A ■ 

where i = 1, 2, 3 sums the spatial components. 

Bose-Einstein condensation can be described in the 
rotating wave approximation, yielding, after normal 
ordering, the Lagrangian density [20] 

£ ma t = gTd^d^ - (m 2 +^R)fti> - M^!,(30) 



where we define the inverse Compton length 

mc 



m - 



n 



(31) 



The equal-time commutator reduces to 

n° $(r, t), doft (r', t)] + n° [ft (r, t), doftr', *)] 
<5< 3 >(r-r') 



ihc- 



-93 



which reads when written out in eigenfunctions 



'-93 



(32) 



(33) 



Multiplying this with yj — gs(r)ip* (r)d 3 r leads to the 
orthonormality (2i?;/7i 2 c 2 ) J d 3 r n°-^—g3ip*tpj = Sij- 
The Hamiltonian density becomes 



(34) 



H = d v ftd ij) + d'ftditp + (m + £R)ftip 

The non-relativistic, flat space Gross-Pitaevskii 
equation reads [20] 



ina t * = -|-V 2 *o + ff|*ol*o, 
2m 



(35) 



with g = A-Kh 2 a s /m modeling the two particle interac- 
tion by the scattering length a s . For hydrogen in flat 
space one has [21] 



a s = 0.32 ao singlet state, 



a s = 1. 34ao triplet. 



We shall continue with the singlet value. Our first task 
is to connect A to g. The relativistic form of (35) is 



rf v dud v ^ + m 2 i/> + ^^|V>o| 2 V>o = 0, 
2nc 



(36) 



where rf v — diag(l, — 1, — 1, —1) is the Minkowski met- 
ric and tpo = Ticfy / V2Norru? ', while the coupling, 
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A = 



8m 2 eg 



= 32n 



a a mc 

h 



= 327ra s m = 8.10-10 6 , (37) 
(38) 



is dimensionless. Eq. (36) derives from the groundstate 
of a field (23) with iVo particles in the groundstate of 
a field theory (30) with metric g^ v = rf v and £ = 0. 
Going to curved space is achieved by taking a gen- 
eral g^ v , which yields the bare Lagrangian density, still 
having £ = and A from (37). Renormalization argu- 
ments will add the £R term, and we shall choose £ such 
that the Bose-Einstein condensed problem has a self- 
consistent solution. Thus, it is assumed that £ is set 
self-consistently, and it may consequently depend on 
BH parameters such as the mass. 

3.1. Self- consistent field theory 

The material energy- momentum tensor (T m )„ = 
(p m ,—p™,—p™,—p™) is derived from the quantum 
field theory. We do not include the effect of the £J?- 
term in it. The energy density reads 

_ {dtftdtfy {dr$drj>) . {d e ftdei>) 

pm — ^77 I 77 h 



c 2 U 



V ' W 2 
W 1 sin inc 



The pressures (p™,p™, p™) have this shape with signa- 
ture (+ + ), (+- + ), and (+-- + ), 

respectively. Spherical symmetry will imply that p™ = 
p™ = p™. For a uniform groundstate p m is isotropic, 



/ x E o it 2, , / mr 1t2 , A|*nl 



(40) 



They consist of a vacuum part p — — p = const, and a 
stiff part p = +p ~ 1/U, the types studied in [6] and in 
[15]. In the non-relativistic (Eo = mc 2 ) and flat space 
(U = 1) limit, they reduce for A = to the expected 
results p m — mc 2 ^^ an d Pm = 0. 

The energy momentum tensor in the Einstein equa- 
tions has two further terms. Because of the ^R-term in 
(30), the Einstein equations embody a direct backreac- 
tion (B) of matter on curvature, G^ v = 8-irGc' 4 (T m v + 
t k + T bT + r iT) - BG^, with 



B = (^V)) = B + B e 



c 

Bq = — 7-t,N\il> \ 2 



(41) 



If B depends on r, the derivation of the Einstein equa- 
tions brings from the term £R(ffiip) also derivatives 
of B, which induces another new term T£" . Let us in 
general define the dimensionless density and pressures 



P = — 7-M p, p r = — —M Pi . 



(42) 



The elements of (T B )i = diag(p s , -pf , -pf, -pf ) 
can then be expressed as 

^2 / B" 2rB' B'U' \ 
p B = M [— + 



(43) 



V W 2 2UV J ' 

B 2 / B'U' 2B'W 

— —M 1 

\2UV + VW 

—2 ( B" 2rB' B'W 



-B 
P± : 



-M 



V W 2 



VW 



where the explicit factors r arise from inserting the 
harmonic constraint (11). 

In terms of the function Bo, eq. (39) and the pres- 
sures read 



p m = 



Pr 



P± = 



E Bo 


+ 


R' 2 
° + 


m 2 


U 


WB 


16tt£ 


E B 
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Bo 


m 2 
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WB 


16tt£ 


ElBo 




R' 2 


rn 2 


U 




WBo 


16tt£ 



where 
7? - E ° 



(B + ^B 2 ) + p e , (44) 

(B + ±B 2 )+p*, (45) 

(Bo + ±B 2 )+pl, (46) 
(47) 



(48) 



and where p e , p% and p\ arise from the excited states. 
They are derived in (168). Together with B e , they rep- 
resent all fluctuations of the problem about the ground 
state at temperature T > 0; in fact, even at T = 
these fluctuations exist and they already bring impor- 
tant effects, see [20] in general and the discussion be- 
low for our BH. Furthermore, the Bq ~ \tpo\ term is of 
relative order 



A 



Amp 



327rm 2 £ ' 

Eq. (8) involves a total energy momentum tensor 



(49) 



-'tot 



"1 ' A ' bi "T" _ B — rpfJLV _|_ rpfJiV _|_ rj-tflV 



l + B 



(50) 
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_„ tot 



with elements (T to t)£ = diag(p to t, — p r ot , 
As mentioned, of the last two terms in (50) only 
the pu/U contributions are relevant for us. = 
diag(p, — p r , —p±, —p±) has dimensionless elements 



+ 



M \ B" 2rB' B'U' 
l + B l~V + W 2 ~ 2 

m 2 m 2 - 2 2B ( 2 B 
+ 2F B ° + 4f ABo + 4£t/ 



8£B V 



M 



Pr = 



M 

l + B 

■2 



B't/' 2B'W 



B 



/ 2 



2C/V 



m " R m2 X R 2 

2^° " 1T^' 



8£B V 
.2 



(51) 







2E 2 B - nijtjB p; 



Px = 



M 

l + B 
■2 



B^ 



2rB' B'W 
+ 



B 



VP 



m- m 2 - 2 2EqB, 
—Bo - — AB + 4 ^ 



PI 
M 2 



The Ricci scalar follows from R — p^-j2U + R„ 
where R m = StyGT/c 4 reads 



R„ 



ai+B) 



3£B" 6£rB' B' 



V 



+ 2rnB a + rn\Bl 



W 2 4B V 

2E 2 B - nUB jf e 
2U + m 2 



(52) 



The last term is proportional to T e — p e — p% — 2p e ± , 
the trace of the energy momentum tensor from excited 
states. 

The Gross-Pitaevskii equation, 
2 

^Vo + ^Vo + if^o = (m 2 + £R)ip + ^^o,(53) 
may be expressed in terms of Bo ~ ipo and B = Bo + 



-6£B„( — + -pps-) - (1 + B)(— + -^5-) 
+ (1 tf e A 5 ° + 2m 2 B [(l + AB )(1 + B e ) - Bo] 



+ - 



2b; 



2b v 

- Atgjg D 2 PbiC RR 2gT e Bp 

"77 ^0 Ft - t>oD e = — = z — . 

U U M 2 



(54) 



We can now first verify that the total energy momen- 
tum tensor is conserved due to the harmonic condition 
(11). The terms B 2 /U are singular because we shall 
consider (7(0) = and Bo(0) > 0; they drop out from 
(51,51,52) and (52,54) for E = p, hi ^jt/2, that is, 



t-,2 1 , t 2 2 2 SuGTl 

E = c H hi = £p bi . 



(55) 



We shall verify later that B e (0) = 0, so the B B e /U 
term is indeed less singular. We may decompose f as 



Co 



(56) 



1+Ce 

We shall later see that £ e is small. The leading term is 

Co s |^ 2 m 2 = ^L_M! = 3.288 • 10 55 M| (57) 
3 imp 

In the regime B, B 3> l/£, Eq. (54) simplifies, 
-2 



Bq 2rB' m 
— ~ W 2 + 

B'J 2rB' m 2 



+ 



(58) 



v w 2 ^ (1 + XBo)Be + ^ + ^' 



In terms of V>o the Gross-Pitaevskii equation (58) reads 



Vo' 2r^ B 2 



V W 2 4B 2 V 



■00 + 



2 [1-(1-A)B ] 



6£B 



tpo 



2Br 



■ipo- 



(59) 



Either (58) or (59) determines Bo or tpo, once B e and 
T e are known. To leading order, they can be omitted. It 
is simpler, however, to work with B — Bo + B e rather 
than with B itself, which satisfies 



B" 2rB' m 2 



3£ 



V W 2 



(l-B + XB) 



2rn 2 M 



35 

It also holds that 



m 2 +£ > Rm = — Xm 2 B 



[2 — A(l — B + B e )]B e 



+ 



: — Am B 



Sb 


(60) 


~2M 2 ' 




MbiBe 


2T e 


3U 


3 ' 


rB' 


(61) 


+ BW 2 


STm 





452 V 4BM 2 6BM 2 ' 

where the first relation is exact and the second holds 
for B > l/£. 



4. Exact solution in the interior 

Since A will turn out to be very small, we first take it 
zero. We also neglect the excitation terms B e and T e . 
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Solving Eq. (60) for A = and a constant B, we find a 
relation and, with (55) as consequence, 



5=1, iV|Vo 2 | 



Mbi 



go _ Pbi_ 

8tt£ ~ £ " 8ttV2?4 



(62) 



The value B = 1 expresses a 100% direct backreaction 
of matter on the metric. 

Instead of searching a finite U, as for boson stars, 
[9,10] we assume a very small U with U(0) = 0, coded 
by a parameter it, of the form 



U = \u^ 2 hi W 2 



(63) 



It is custumary to introduce the mass function M(r) 
and M(r) = GM(r)/c 2 , defined by 



V : 



(64) 



1-2M/W 

The 00 and 11 Einstein equations then take the form 

rj-ji 4-xG l 2 
M = -^VK W ptot, 

^-M) U --^ = ^W'^. (65) 

We may combine (10) and (51) together with B = 1 
and the Ansatz (63), to obtain 



(ptot.Pr *) 



1 , m 2 W 2 



4nGW 2 \ 4u 8£ 
In RTG there is a solution of (65) and (66), 

t-j w m 2 w 3 



(66) 



(67) 



It would not exist within GTR, as is seen by taking 
Mbi — » 0, u — > oo first. For the Schwarzschild black hole 
the horizon occurs when M = M for W = 2M. Con- 
cerning the outside metric, we will be close to that sit- 
uation. This implies that a mass M requires the large, 
non-constant coupling £ = £o£i with 



£ = ^1 = 3.288 • 10- 



( M 



3m 4 p ' I 10 9 M G 



(68) 



and £i = 1 up to possible small corrections to be dis- 
cussed further on. Interestingly, £o is a measure for the 

2 

area of the black hole (4irM ) expressed in units of the 
square of the hydrogen Compton length 1/m = h/mc, 
so it is a measure of the surface entropy. Notice that in 



the Bekenstein-Hawking entropy the area is expressed 
in squares of the Planck length. In this work we shall 
consider the T = situation, so that the entropy is 
exactly zero. 
We now get 

— 1 1 

i£ = —pbiGMm = —p, hi mc 2 . (69) 
V 3 V 3 

and coupling 
A = 



\m 2 P 3a s h 2 c 4 _ A -ice ln -i3 1 



327rm 2 £ 2G 3 m 3 M 2 4 ' 166 ' 10 M | • 



Though the latter is small, the product A£ is large, 



A = 6A£ = ^ = 8.19-10 43 . 
Gm 



(70) 



Let us introduce the 'Riemann' variables x and y by 
W 



vv r 
2M y V 



(71) 



so that (7 = 2pl i x 2 . With (63), (64) and M = \M(x + 
x 3 ) from (67), the harmonic constraint (11) brings 



8rx'- 



2x' 2x" 2xx' | 4x' _ 
x x' 1 — X 2 X x 2 (l — x 2 )' 

Going to the inverse function r(x) makes it linear, 
x 2 (l - x 2 )r xx + x(3 - 4x 2 )r x = 4r. (72) 
In terms of the variable y this transforms into 

x4 ryy ~ ^x 2 yr y = 4r. (73) 
The solution is then remarkably simple, 

V5/2 



r= n (l + JL) l^y 



(74) 



(The second independent solution with y/Z — > —a/5 or 
y — > — y is singular atr = 0,y = l.)lt will hold that 
ri ~ M. This determines the metric functions W' and 
V, 



W' = ^x 2 -^y(l + y)^, 



(75) 



r 2 

Putting these results together, it now follows that 
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p = 



3c 4 



Pr = P± = V = ~- 



3c 4 
64ttGM 2 



(76) 



So in the interior we reproduce the vacuum equation 
of state p — ~p = const., that is, p — —p = 3/8. 

To understand the structure of the problem, we again 
take A = 0. Then for any A there is the solution 



B(x) = 1 + Ay = 1 + Ayjl -x 2 . 



(77) 



Surprisingly, their A-dependence factors out, keeping 
a vacuum equation of state p = — p = 3/8, so (77) is 
an exact, non-uniform solution of the same metric. It 
can be verified that ipo ~ y/1 + Ay solves the Gross- 
Pitaevskii equation (60) at A = and without excited 
states terms, as it should. 

4.1. Normalization 



i 

N = No + N e = 2^3— / Ay [B (y) + B e (y)]. (81) 
m J 



It allows to relate the BH energy to the groundstate 
occupation, 



Mc = vNomc 2 , v 



(2 + A)VT 



(82) 



Clearly, the energy Mc 2 of the BH can be at best 29% 
of the rest energy Nome of the constituent hydrogen 
atoms. In the BH formation, the major part of the 
energy, {l — v)Nomc 2 minus the potential energy, has to 
be radiated out. This may explain the large luminosity 
of quasars. 

The leading part of the energy of the quantum field, 



To normalize i[>o, we need the 3d volume element in 
the future time direction, dS M = drd#d0n M ^/— g A = 
SgdV, set by the timelike unit vector n M = S^/VU 
and g3 = ~VW 4 sin 2 6>. This results in dV = 
drdQy/V/U W 2 . It then holds that ^/=gd 4 x = 
Udx°dV. 

The general inner product [19] (ipi,tp2) = 
(—i/hc) J dS M (ipidftipZ — <9 M y>ii/>|) defines the or- 
thonormality, already noticed below (33) , 



I 



dVipitpj = 5i : 



(78) 



together with (ipi,tp*) — and (ip*,ipj) = —Sij. With 
the volume element 

dV = dydn8M 2 /Mbi 

it yields a volume V = 327rM 2 //Xbi 
and 



(79) 

47r(2M) 3 /Mbi 



\/3/i b iC 



.6 



l*Sl - |§ivc|^l = l67rG . imM2 



(l + Ay). 



(80) 



4.2. Properties of the solution 



equals 



E {0) =- - 
* 4/i bi G 4 m 3 M 2 



(84) 



With \/p bi = 17.5/M| it is of the order of the total 
energy Mc 2 , and it may exceed it. To obtain the to- 
tal energy, the gravitational energy density has to be 
taken into account, which is partly positive and partly 
negative. Due to a sum rule the total energy is always 
Mc 2 [2]. 

At the origin the solution exhibits the powerlaw sin- 
gularities 

U = U ir ^, V = ^lW 2 r~<»- 2 , W = WiA" 1 " ,(85) 

where 7 M = ^(v^ + 1) is the golden mean. But if we 
take W as the coordinate, we have in the interior the 
shape 



The horizon is located at r = M , x = 1, y = 0. We 
may integrate (80) over the BH, which yields iVo- Alter- 
natively, from the definition (41) of B, we may consider 
J dV B, making use of the normalization (ipo, ipo) = 1. 
Either way, this yields 



ds 2 = ^ /j, 2 iW 2 c 2 dt 2 



2dW 2 



1 - IV 2 /4M 2 



W 2 dD. 2 , (86) 



which is regular at its origin, with the term 2dW cod- 
ing the above powerlaw singularities in r. 
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5. General solution near the horizon 



In the region B>l/^we have the simplifications 



Near the horizon the exact solution will be deformed. 
In general we may code the functions U(r) and V(r) 
in new functions u(r) and v(r), 



U = 2fi bi x u, V — 



8M x' 2 

2 1 

y 



(87) 



and W — 2Mx. With (11) we then have for arbitrary / 

frr 2rf r _ X f yy \_a , X ( U V _ V V \] fv 

V + W 2 - 8M 2 v 1 V 2 l - U ,\mV 

since the W" terms cancel, as they should, because 
one can also start with W as variable instead of r. The 
leading shape of the Gross Pitaevskii Eq. (60) can be 
written as function of y, (we neglect A, as it will be 
much smaller than B e and T e ) 

2 2 

X -B vv - y-B y + |-(^ - ^)B y + B - 1 = S% 

S% = 2(1 + &B e - ^B e - |f e - UB - 1). (88) 

In this equation B — B + B e consists of contributions 
of both the groundstate and excited states, while B e 
and T e involve the latter only. While v = 1 in the 
BH interior, it grows as ~ /v beyond the horizon, 
possibly enhancing the effect of the flucutations. In 
principle this equation may therefore describe a decay 
of Bo to zero, embedded in excited states that decay 
slower. 

The Einstein equations read in terms of the new 
functions u and v 



\i2- y -) + 1(2-2*; 
x z v v v 

2 ^- 2 



(89) 



X z X Z V UV a+tt 



Expressing the shapes (51) in y, we have 
__x 2 B yy -(3y+^)B y 6S +3AS 2 



8v(l + B) 



- _ (3g-^)g» 65 + 3A5 2 . 



86(1 + S) l + B' 
(90) 



P± = ■ 



8v(l + B) 8(1 + 5) 



8u(l + B) 



86(1 + 5) 



Pi 
l + B' 



-_ (v- E £ L )Bv , 4 + 25 + A(45 + 35 2 ) . 



P± 



8v(l + B) 8(1 + 5) 

-ySy 4 + 2B + A(45 + 35 2 ) 



8v(l + B) 86(1 + 5) 



+ 

with the following source terms at A = 
6(2 + B) + 4p e - 3(1 + 6)5 e + 2S| 



Pe = 



4(1 + 5) 



-» - _* 6B J. ^% + 3(1 + £ e )Be 



8(1 + 5) 4(1 + 5) 



(91) 



The equation (72) for the Minkowski coordinate now 
reads 



2 2 , 

x y t X x T" £ 



4y 2 _ i + £?L(^ _ ^) 

y 2 v it u ' 



r x = Avr. (92) 



6. Excitations 

Let us reformulate our theory on a new basis. We go 
to a new coordinate z, 



W 1 z 

x — — =r- , y = — tanh — = , 

2M cosh( Z /^2) n/2 



(93) 



so that z — — oo at r = and z = at r = M, a 
dimensionless time s and a scaled energy, 



V2 rac^bi 
The line element then becomes 

ds 2 = 4MV(uds 2 - vdz 2 - dfi 2 ) 



(94) 



(95) 



which is Minkovskian in the exact solution where u = 
v = 1. It corresponds to volume elements 

d 4 r-/=<j = d 4 f v/w 16M V, d 4 f = dsdzdft (96) 

„, , ,„8M 2 /F , Jn 4M 2 i 2 /2w . n _, 

dV^dydfi x -=dzdQ \ — (97) 

fJ-bi V u /i bi V u 

In terms of the new field 

8^^^-a^ = \ X (98) 
the innerproduct (8%l>i,5ipj) just becomes 
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( Xi ,Xj) = {Ei + E 5 ) J Az&Sl-xlxi = Sij. (99) 
For the kinetic term it holds that 

d 4 r y /=gd' 1 6ftd li 6ii> = d 4 fvhc ( MlM (iqo) 



1/2 2 \ 



where Xi = ( V / U ) 1 ^ 4 X./ X an d L = — i[de, (1/ sm9)d^\ 
is the angular momentum operator in units of h. Like- 
wise, the — (m 2 + ^Rjd'tp' Sip term becomes 



w 16M a; 4 Vu 



(101) 



with potential 



1 s 2 2:r 2 T e 



3B 



In vacuum, 2: — > + , i? m — > and Vbi — > 6£, which acts 
as an infinite barrier. More precisely, it leads to a decay 
at a scale of the Compton wavelength h/mc = 1/fn. 
Finally the interaction term brings quadratic terms 
from the expansion ip = aoV'o + Sip up to second order. 
Noting the time-dependence ipo = |V'o|e~ lMS with 



this yields in a straightforward manner 



(102) 



=4 — (4x X + e M X +e ^ X ), (103) 

16M 

where Bo is the groundstate contribution to B. Notice 
that (101) and (103) involve the same prefactor A£o = 
A/6£i. Together this brings an action 



S mat = hj dsdzdnv(L 2 + L int ) (104) 
with, after a partial integration, a Lagrangian density 

L2= dAx_l dz ^ dzX + L ^. Lx (105) 



(^ + V + Vi)tfx - \Vr{e 2 ^x 2 + e" 2 ^ 2 ) 



' y u ^ "' ' 1 ' ' ~ - 
which involves the potentials 

,_4 D 2 

vb( z ; 



— -(l-A--)x - — 



l-2a; 2 y V i/ 
2v 2\/2 « « 



(106) 



16t) 



4m" _ 4t/\ _ V _ t/ 3t/ 5v\ 
W M u « u V ' 



K ' cosh 2 (z/V2) 6 

The equation of motion for the x field is 

_9|X = + ( i + L 2 + V() + y i)± + Vie -2iHs x l 

U V u 

We have to perform a Bogoliubov transformation in 
our case with non-constant "potential" £R, a situation 
similar to the case of non-constant external potential, 
discussed in e. g. [20]. With i = (n,£,m) we set 5 

X (z, 9,4) =^ [ u t {z)Y em {6, fle-'V+^'bi (107) 

i 

+ v;(z)YU0,<t>)e- i{ *- Wi)s bt] , 

and the conjugate momentum is (v/u)d s X- With f = 
(z, 9, (f>) the commutation relation (32) becomes 

+ "OMfW) + «?(*)**(*')] (108) 

i 

- - - Wi)[ui(f )«*(?') + m* (?>;(?')] = S(r - r). 



where w»(f) = Ui(z)Ye m {9 , <f>). We may then expect the 
orthonormality property 



/ 



dz ^ [u1(z)u 3 (z) - Vi{z)Vj( Z )] = 



(109) 



The eigenmodes satisfy the Popov equations 



2p,LJi u'l 

Mi = h 

M V 



V + £(£+!) 



+ V 1 (u l + Vi), 



2pWi v" 
— Vi = h 

U V 



Vq +£(£+1) i - 



+V 1 (u i + v i ), 



(110) 



where derivatives are with respect to z. From these 
equations one can show that 



(wi — oj*) / dz-(uiU* — ViV*) = 



(111) 



5 The complex Bogoliubov functions Uj and i;, should not 
be mistaken for the positive metric functions u and v. 
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and with real w; this confirms the inner product (109). 

In the interior one has u = v — 1. At uii = I = 
and T e — > the above equations then allow the exact 
solution Mo = —vo = z-y/l + At/, which corresponds to 
±•00 on this basis, so the solvability stems from the 
one of the Gross-Pitaevskii equation. This situation is 
related a gauge transformation that changes the phase 
of the groundstate ipo [20] . 

For excited states ~ £ ~ V ~ A are large, so 
both Vb and u> 2 can be neglected. Then lj ~ A/L In the 
regime 2: <g; — 1 one has V <C 1, so one expects 



ik n z 



-Ve 1 



(112) 



We can deal with the boundary conditions of these 
excited states as with plane waves, e. g. by requiring 
that e lk " z = 1 at some large z = — L and take L — > 00 
at the end. The normalization constant is then 

(113) 



Another regularization is to assume that k = k' — ik" 
has a small imaginary part; this will keep all integrals 
starting at —00 finite. Then c n ~ \fkF . Below we shall 
employ sine-modes and impose the hard wall boundary 
condition at 2 = —L. 

6.1. Between the center and the peak of the potential 

In the typical case where A > the potential V = 
A(l — y 2 )(l + Ay) has a maximum at 



r- , y/l + 3A 2 - 1 
z c = -V2arctanhj/ C , y c = — , (114) 

which goes to zero for A — ► 0, but remains finite for 
A — > 00. The region — 00 < 2 < z c , which covers the 
whole interior when A — > 0, is considered first. 

Since V = Xx 2 B is large, the excited states can be 
analyzed with the WKB method. The function 



S = Ui + Vi, 

satisfies with 



the equation 



(115) 



(116) 



k 4 s = s lv - 2£ 2 s" + ts + ^ 2 2\/s - 2{Vs)". (117) 
We make the Ansatz 



s(z) = const, e 



iS(z)-T(z) 



(118) 



where S = 0(A 1 ^ 2 ), r = O(A ) and higher order cor- 
rections may be neglected. The role of L will be dis- 
cussed below. We have at leading order 



S' 4 + 2(V + P)S' 2 + t + 2V£ 2 = k 4 . 



with the solution 



S' 2 = yjk* + V 2 -V -I 2 



(119) 



(120) 



At a given location z the solution is of plane wave type 
when S' is real, which occurs provided I is limited, 



t 2 < l\(z) = ^k 4 + V 2 {z) - V(z) (121) 

We shall not need S itself. At next order we find 

4(S"' 3 - S'V - S'£ 2 )t' = 2(S' :i - S'V - S'l 2 )'. (122) 

with solution fixed to r = for z — > —00 (V — » 0), 

1, k 4 + V 2 1, Vfc 4 + V 4 — V — I 2 
r= 4 ln ^^ + 4 ln W~I 2 ' (123) 

To get Ui and Vi to leading order is now easy. Since 
(I 2 + 2V)s - s" _ V + Vk 4 + V 2 



Ui — Vi = 



k 2 k 2 
and imposing S(Zi) = 0, we end up with 

1 k 2 ± (V + Vk 4 + V 2 ) 



-s, (124) 



(Ui,Vi) = 



fiL 2k(k 4 + V 2 ) 1 / 4 

k 2 -£ 2 \ 1/4 • • 
x I , smS. 

Vk 4 + V 2 - V - I 2 



(125) 



For z — » —00 (V — > 0) one has indeed m — » 
sin kz/^/jlL, w» — » 0. If we impose a hard wall bound- 
ary condition Ui — 0, Vi = at 2 = — L, the normaliza- 
tion (109) is satisfied, since it is determined by values 
2 <JC — 1. As costumary for plane wave problems, L 
will be taken to infinity at the end. 

For £ > l+(z), the action S becomes imaginary, ex- 
pressing a damping of the wave that has to penetrate 
the potential barrier to reach this position z. These 
states lead to negligible corrections. 

We can calculate the Hamiltonian 
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H = const. + V" Efb\bi, 



(126) 



with 



= J cb[(Kl 2 + l^r+^ 2 +? + vo(kl 2 + M 2 ) 

+ V( Ui v* +u* Vi )]. (127) 
With K| 2 + \vl\ 2 = S" 2 (M 2 + l^l 2 ), this becomes 

El = J dz[^ + ^k^ + V 2 )(\u,\ 2 + \v t \ 2 ) 

+ V(uiV* +u*Vi)]. (128) 

As in the normalization, these integrals are dominated 
by large negative z-values, where V — > 0, so that we 
simply get 



Ei = VC + fc 2 ~ A + 



2/i' 



(129) 



It is degenerate (independent of £, m) as in the quan- 
tum Hall effect, though here there are no spectral gaps. 
At a given location z the wavefunctions m,Vi are os- 
cillating (not damped) provided I < £+{z). 

6.1.1. Contribution to the fraction of excited states 

We can now consider B e , the excited states contri- 
bution to the direct back reaction B = B + B e , that 
arises due to the interaction term, even at T = 0. We 
start from the definition 



and express this as 



(130) 



3m 2 p x 2 ^/u 



(X f x) = eB e (z), 



47rm 2 A 3 / 2 
3m 2 , V? ' ~° A 3 / 2 a; 2 0Z^« 



Be = 



2^ V- 



(131) 



In the exactly solvable case the sine modes, that exist 
for z < z c , i. e., to the left of the peak of V, behave as 
sin fez for z — » — oo. The hard wall boundary condition 
Ui = 0, Vi = at z — L brings the quantization k n — 
niv/L, n — 1, 2, • • •. Replacing sin 2 S — » |, this yields 
at a given position z, 



B e = J^-- I dk 

X i/2 X 2 n 



d£ 2 \vi\ 2 



oo 

/ 



dfc (K 2 + V - k 2 )' 



ksjK 2 - V 



(132) 



+ (k 2 + V - A- 2 )arcsinh 



K 2 - V 



k 2 + V -K 2 

where V = V(z) and K 4 = k 4 + V 2 . The integral gives 

- = °- 236792 y3/2 = re ^^3/2 (133) 



so the result neatly vanishes at the origin (a; — » 0), 



Be = 0.757 B 3/ V 



M 



(134) 



1O 9 M ' 

According to (131) is has the characteristic strength 



47rm 2 A 3/2 _ 247ra 3 / 2 c 2 
3m^ ~ G 3 / 2 m 1 / 2 M 



0.382 



10 9 M, 



o 



M 



(135) 



In standard BEC, the correction to the groundstate 
energy is of relative order \Jna% [20] . Using (82) we find 
in the homogeneous case A = 



\Zn H a 3 s = 0.0136 



1Q 9 M 
M 



(136) 



so our variable e is of the same order of magnitude, 
confirming the expectation of the introduction that the 
relevant physical parameter is na 3 , which is small for 



6.2. Between the peak of the potential and the horizon 

We take V = oo beyond the horizon. Then it has 
quasi-bound states in the region z c < z < 0; they are 
not true bound states because V drops to zero for 2 << 
z c - But since the energy barrier ~ A will be very large, 
the tunneling into the interior will be extremely small, 
and we shall neglect it. 

For U — V = 1 we can now copy previous solution 
(120), (125). At given value of £, a real valued S" starts 
at 2 = Zi set by 



V{zi) 



2£ 2 



(137) 



The surface state "lives" in the interval z% < 2 < 0. 
The smallest 2, arises when z, = z c and it has £ 2 = 
Vfc 4 + V 2 — V c , while Zi — > for the maximum £\ — 
sfk 4 + A 2 — A. The solution may now be written as 
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1 k 2 ± (Vfc 4 + V 2 + V) 



fiL t 2k(k 4 + V 2 ) 1 / 4 
sin S(z) 
{Vk 4 + V 2 -V-£ 2 ) 1 / 4 ' 



(138) 



which also depends on z through V(z). Since sin 2 S(z) 
oscillates fast, it can be replaced by ^, so the normal- 
ization (109) is achieved by 



dz ( Vk 4 + V 2 + V) 



\(k 4 + v 2 )(Vk 4 + v 2 - v - e 2 )} 1 /' 2 



(139) 



Recalling that S(zi) = 0, the hard wall boundary con- 
dition m — Vi = s — at z — can be fulfilled pro- 
vided the phase S(0) is an integer n times n. At given 
£ this defines the eigenvalue k n . 

6.2.1. Contribution to the fraction of excited states 

It is instructive to investigate whether these states 
cause divergent effects for states localized close to the 
horizon, those with \zi\ -C 1. At fixed k and £, we have 



S = k — X + XAz -t = XA(z - zi) 
with 



fc = (fc 4 + P) 1/4 , 
and 



A : 



A(k 2 - A) 
V2k 2 ' 



Zi = 



X + t 2 -k 2 



\A 



The maximal £ at a given z is 



l\{z) = 
This brings 

(k 2 + \h 



A + XAz. 



Li = 



(k 2 +X)y / k 2 -X-P 



k 2 \/XA 
From S(zi) = we get 



k 2 XA 



S(z) = ^VXA(z-z t ) 3/2 



implying 

So we may set 



2n3/2 



(140) 



(141) 



(142) 



(143) 



(144) 



(145) 



(146) 



S- 1 ^-— Vk^X-7 2 . (147) 

dfc 7T dk % XAk 2 

We can now calculate at given small z 



]>> 2 ( 2 ) = J dk J 



dt-vfiz) 



(148) 



di 



(A + P-fc 2 ) 2 



16njlk 2 (X + k 2 )^/£ 2 + (z)-e 2 



where we replaced sin 2 S(z) by \ . Since the singularity 
at £+(z) can be integrated, these states brings no spe- 
cially large contribution near the horizon y = 0. The 
outcome is of order X 3 ^ 2 f^ 1 ^ 2 , as for the modes near 
the origin, so both type of modes bring comparable ex- 
citations, B e ~ e, as one would expect. 



6.3. States close to the horizon 

States localized close to the horizon can be studied 
analytically, we define 



k 4 =£ 4 + 2£ 2 X(l-^), Zi = e -J^ 



z = Cz 

and supposing s(z) = f(zt — z) 



(149) 



(150) 



C 4 f v - 2(£ 2 + X)C 2 f" + 2£ 2 X^—^-f = (151) 

C V 2 



so to leading order 



The solution is 

fi(z)=Ai{zi-z) (153) 

which oscillates for z% < z < and decays for z < z t . 
The combination 



1/3 



, , C Ai 2 (zi - z) TYC'Ai 2 (zi - z) 



2 v ^/^dyAi 2 (y) ' 2^Ti 
is normalized to J dz f 2 = 1, so this results in 



(154) 
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, 2 2 (w + V±k 2 ) 2 2 
The zeros of Ai(^) occur at 



/37rn\ 2 / 3 



V 2 J 
so this yields 



dn = —\J—Zi&Zi 

TV 



n - _±_|^(")|3/2 



(155) 



(156) 



(157) 



Putting things together yields 
o 



B e (z) = 



-Ai z 



2(2\A 2 ) 

(X + 1 - x) 2 
2X(X + 1) 

where 

X = yjx 2 + 1, Z : 



— / d* / 

4 2 )ve y y 



da; 



X-lX 1 / 3 



A 1/?3 z 

21/6^1/3- 



(158) 



(159) 



In the 2, integral only the small values are reliable and 
the result is not of order unity but of order A" 1 / 6 = 4.7- 
1CT 8 . Still, the result vanishes exactly at z = before 
we pass from a sum to the integral. The derivative is 
well defined, however, 



~ x l/6 



24/3^2/3 



OO 



2/3 



(160) 



(X + l-x) 2 2 
2X(X + 1) 

At 2 = it takes the value 

-, 0.157682 A 1 / 6 
B ^°) Jyl > 

while for z — ► — oo it decays as 

5 , 93V 3 r ( j) A ^ 

82V40F ~ x \z\V 6 

where we used 

oo 

/j 5/2 a • 2 / \ ^ //3 r(^) 
dyy Ai (y)= = 



0.03284 A 1/2 

\\z\y fi 



0.471436. 



(161) 



, (162) 



(163) 



6.4. Fluctuation energy of the matter field 

Let us calculate the energy density and pressures of 
the quantum field. We need the following contributions 



Pi = 



P2 



P 3 = 



he 



° 2U 16MVu 

(8^L 2 Sj>) 
W 2 



16M x 4 v 



he 



16M a; 4 



(164) 



AM x 2 



3hc£ 
Ud^x 2 



Pl = 



4ftc 

3\Bhc£ IA x X ^ a -2 lf ,s^2 + ^^2 



16M x 2 
They determine 



2 >- 



Pe = Pi + P2 + P-A + Pa + P5, 

e e . e e e e 

Pr = Pi + P2 - Pz - Pa - Ps, 

e e e e e 

P± = Pi ~ P2 ~ Pa - Ps, 
which implies that 

T e = -2p\ + 2p\ + 2pt + Ap\ + 4p|. 



(165) 



(166) 



The leading terms at T — are indeed also of order e, 



Pi = 



8a- 2 u 



pt = ^Be, 



(167) 



the other terms are smaller by a factor A at least. So 



Pe 



1 



4uar 



+ 3)B e 



(168) 



This now implies that the sources of the GP equation, 
(88), and of the Einstein equations, (90) and (91), are 

S B = 0, p e - jB e =p% + jB e = B e (169) 



6.5. Reaction of the metric on the fluctuations 

In order to investigate whether the matching of in- 
terior and exterior solutions can be achieved near the 
horizon, we investigate the reaction of the metric to 
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fluctuations of the exact solution caused by the terms 
p e , p% and p e ± . We set 

B = l+Bi, w = l + ui, v = l + vi, £1 = l-p&Q) 
where Bi, ui, Mi and £2 are of order e. 

6.5.1. TTie solvable case A—0 

The B-equation becomes to linear order in e 

^B'i -yBi+B^—n. (171) 

In this section, derivatives are with respect to y. The 
homogeneous solutions are 



BP = y, 



_L _ 1 

2x 2 2 



B W = ^ _ r + r*in 



3j/ 

4 



1 + 1/ 



and they have a Wronskian 

W = B 1 1) B 1 2) '-B 1 1), B 1 2) = ^. 



The solution for Bi therefore reads 



(172) 



(173) 



Si = 613/ + B 



sf^y) / Ay 



?\y) J 



dy 



8ST m B\ 



(2) 



3x 2 W 
88f m B { ^ 



3x 2 W 



+ ^B 1 2) ( 2/ ) / dxx 3 ST, 



biy 



(174) 



^ / A X X -Sf m Bf\ 



with 61 an integration constant. The last expression 
exhibts regularity at the origin x = 0, even when 5T m 
has al/i singularity, as we discussed above, 
tii may be solved from the }- Einstein equation, 



Vl 



ui — x yu 1 
3x~ 2 



UpT + HBi-yB'!) 



after which mi satisfies 

2 11 . 1 . 4 
x Mi - 4yuj + — wi = 
x 2 

with source term 



-4pi + 14p 4 - 75- (Pi - p 4 ) 



(175) 



(176) 



(177) 



The homogeneous solutions P/x and Q/x involve the 
associated Legendre functions P = P^, Q = Qi^, 



The solution then reads 



Ml 



y v 

Ps P { , Qs 



dy xW ~ x / dV xW 



- I 

X 



Q± _ Q f dx £± 

yW x J yW 




with the Wronskian 



W = PQ' - QP' 



8i 



V3i 



(179) 



(180) 



It is imaginary, so mi is real. (179) is regular at x = 
0, so no homogeneous solutions can be added. It now 
follows that vi(y) diverges as C-2/y 2 near the horizon 
y = 0. We get 

c-2 = ^ - |[jJi(0) - j5 4 (0)] + iSi(0) (181) 



with 



Ml(0) 



dw(Q - B)xs 



(182) 



(183) 



For A = one has 

Pi = — , p4 = 6p a a; 
a; 

with p a = 2.80959- 10" 15 M| a positive amplitude. This 
gives mi(0) = 13.6518 p a , while Bi(0) = -496p a /45. 
Together they yield 



c_ 2 = 5.70616p a 



(184) 



which is positive, and showing an upturn of v in the 
narrow region y ~ near the horizon. 

6.5.2. The general case A > 

When ^4 > 0, the Bi equation gets coupled to the 
ui , vi equations. No explicit solution of the linearized 
problem has been found. Inspection of the equations 
near y — reveals that now a singularity 



+ O(y ), (185) 



vi = —2-+0(y), Si = — ''' ' 

2/ 2 2y 

is allowed. From the above case A = it is to be ex- 
pected that c_2 remains positive, so at A > also Bi is 
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singular. The signs are the ones expected for approach- 
ing the vacuum: B decays while v increases, towards its 
high peak slightly beyond the horizon. In retrospect, 
the induced decay of B also indicates that A > is the 
typical case, rather than the no-hair value A — 0. This 
important fact gives hope that a self-consistent treat- 
ment near the horizon achieves to match the exact so- 
lution in the interior with the deformed Schwarzschild 
metric in the exterior. It remains an open problem to 
consider this behavior in a self-consistent way. 



7. Conclusion 

We have questioned the general wisdom that static 
BHs have all their mass in the center and that its inte- 
rior cannot be described by present theories based on 
General Theory of Relativity. Estimates show that a 
picture of closely packed H atoms naturally applies to 
the supermassive BH's in the center of galaxies, M ~ 
1O 9 M . We therefore attempt to describe them as more 
or less normal objects like stars. 

We present within the Relativistic Theory of Grav- 
itation (RTG), an exact solution for a BH, of which 
the interior is governed by quantum matter in its Bose- 
Einstein condensed phase. Powerlaw singularities oc- 
cur at the origin, that get absorbed in the Riemann de- 
scription of the metric. Elsewhere, the solution is regu- 
lar. The redshift at the horizon is finite, though of the 
order l/£ bi ~ 10 14 . 

This solution is still to be matched with the 
Schwarzschild metric, which near the horizon is de- 
formed in RTG. We have carefully derived the com- 
plete fluctuation spectrum about the groundstate. 
The matching of the inner metric with the outer met- 
ric at the horizon has been considered in a first order 
perturbative approach, which shows an enhancement 
effect near the horizon. The full problem still has to be 
carried out, and it has to be done self-consistently. It 
remains as a task for future to show that this indeed 
leads to a proper decay of matter and behavior of the 
metric near the horizon. 

Our BH is a quantum fluid confined by its own grav- 
itation. It puts forward that a BH is just an intense 
gas cloud, without an event horizon, as was also de- 
duced from the Schild et al. observations [1,3]. In the 
interior, time keeps its standard role. No Planckian 



physics is involved; Hawking radiation is absent and 
Bekenstein-Hawking entropy plays no role. In our zero- 
temperature situation the entropy of the quantum field 
vanishes. Because the Schwarzschild singularity is cut 
off by the bimetric coupling, there is no connection with 
any form of quantum gravity, even though the redshift 
at the horizon is of order 10 14 . 

Schild et al. [1,3] have explained their observations 
in term of a magnetic dipole moment of the black hole, a 
"hair" . Our BH also has one "hair" , the binding energy, 
expressed as -Ebind = Nrnc 2 — Mc 2 = (1 — u)Nmc 2 . 
Here v can take any value below l/2\/3 = 29%. We 
confirmed that the previously derived solution of the 
Gross-Pitaevskii equation with a free parameter, A or 
v(A), shows up also as a zero mode in the fluctuation 
spectrum. As one would expect for a classical theory 
of gravitation, when the quantum matter in the BH 
has reached a certain groundstate, the classical metric 
allows the system still to go to a lower energy state. 
Indeed, the passage of celestial bodie will induce os- 
cillations in the metric and emission of gravitational 
waves, which, upon re-equilibration, increase the bind- 
ing energy, finally up to 100% of the rest energy of 
its constituents, Nmc 2 . In that final state the mass is 
completely balanced by the binding energy, making it 
look like a zero mass object. In its stable state, the BH 
has a fraction strictly-less-than-one-half of the ground 
state energy of the constituents, so the major part of 
the zero-point energy has to be emitted in radiation. 
This property may explain the enormous luminosity of 
quasars. 

It has been assumed that renormalization couples to 
the matter field density to the curvature scalar with 
a strength £. This parameter is chosen appropriately 
in a first step of renormalization of the scalar field 
theory with a large quartic coupling. The value £ = 
2m 2 M 2 /3m 4 P ~ 3 • 1O 55 (M/1O 9 M ) 2 allows an exact 
and explicit solution of the interior metric, any other 
value would be inconsistent. The fluctuation spectrum 
is well defined, and various very large or very small 
numbers finally combine into reasonable prefactors. It 
is noticed that the leading corrections of the matter 
field are of order V naj , as it happens for in Bose- 
Einstein condensation in a box. The physical reason for 
this, diluteness of the gas of H-atoms, was put forward 
in the introduction. For this reason we expect that the 
field theory for supermassive BHs with M 3> 1O 9 M 
can be renormalized perturbatively after the first step 
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that fixes the leading value of £. 

We have set out the lines for studying the fluctua- 
tion spectrum near the horizon. It is left as a task for 
future to show that they indeed fluently connect the 
empty space metric of the exterior (i. e. the deformed 
Schwarzschild metric) with our exact solution for the 
metric in the interior. 

An important question is whether formation of re- 
alistic supermassive BHs brings the matter indeed in 
or near the Bose- Einstein condensed groundstate. This 
would require the study of the finite temperature situ- 
ation. Extension to finite temperatures, not presented 
here, will exhibit a y 3 / 2 fraction of thermal atoms. Also 
the stability of the solution needs to be studied. 

Calculation of the normal mode spectrum may lead 
to predictions that deviate from the ones of GTR; this 
spectrum may be observed in the foreseeable future. 

We failed to apply our approach to GTR, technically 
because it lacks compensation for the 1/U terms, that 
in this solution are truly singular at the origin. If no 
other solution exists for the considered physical situ- 
ation, GTR must be abandoned and replaced by an- 
other theory, RTG being the first candidate. In view of 
its smaller symmetry group, this may have far reach- 
ing consequences for singularities in classical gravita- 
tion - they would probably be regularized - and for 
quantum approaches to gravitation, since the primary 
space-time, Minkowski space-time, needs no quantiza- 
tion. 
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